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Abstract 

In the framework of the minimal cosmological standard model, the ACDM 
model, the Dark Matter density is now known with an error of a few percent; this 
error is expected to shrink even further once PLANCK data are analyzed. Match- 
ing this precision by theoretical calculations implies that at least leading radiative 
corrections to the annihilation cross section of the dark matter particles have to be 
included. Here we compute one kind of large corrections in the context of the min- 
imal supersymmetric extension of the Standard Model: corrections associated with 
two-point function corrections on chargino and neutralino (collectively denoted by 
X) lines. These can be described by effective x~fermion-sfermion and x~X~Higgs 
couplings. We also employ one-loop corrected x masses, using a recently devel- 
oped version of the on-shell renormalization scheme. The resulting correction to 
the predicted Dark Matter density depends strongly on parameter space, but can 
easily reach 3%. 
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1 Introduction 



The presence of Dark Matter is by now well established. In fact, combining data from 
several cosmological observations, including in particular analyses of cosmic microwave 
background anisotropies, allows to determine the universally averaged Dark Matter den- 
sity with an error of about 3% [T], assuming that the standard cosmological "ACDM" 
model, which combines a cosmological constant with essentially non-interacting cold Dark 
Matter, is correct. Note that this simple model describes all cosmological data satisfac- 
torily. This error is expected to shrink to the 1% range once the PLANCK collaboration 
[2] releases their data. 

A particularly attractive class of particle physics candidates for Dark Matter are 
Weakly Interacting Massive Particles (WIMPs). Their relic density can be calculated 
precisely in the framework of standard cosmology [3] , assuming that the reheat temper- 
ature after inflation was higher than about 5% of the WIMP mass [H |5j |6]. Input in this 
calculation is the WIMP annihilation cross section into Standard Model (SM) particles; 
the predicted relic density scales essentially like the inverse of this cross section. A pre- 
cise prediction for the relic density therefore requires an accurate calculation of the total 
WIMP annihilation cross section. 

This can only be accomplished within a specific WIMP model. Since the SM doesn't 
contain a WIMP (nor any other plausible candidate for cold Dark Matter), it needs to 
be extended. The best motivated WIMPs are those that arise in extensions of the SM 
that have independent motivation, not related to Dark Matter. 

Here we focus on supersymmetry, more specifically, on the minimal supersymmetric 
extension of the SM (MSSM) [7j. Superpartners of all SM particles not only stabilize the 
weak hierarchy against radiative corrections, they also allow single step Grand Unifica- 
tion of all gauge interactions. We assume that R— parity is exact, so that the lightest 
superparticle (LSP) is stable, thereby satisfying one of the conditions on a successful 
WIMP candidate. Moreover, we assume that the supersymmetry breaking terms con- 
serve CP. Barring several finetunings, experimental constraints on CP violation, mainly 
from electric dipole moment measurements, require their phases to be small, unless su- 
perpartners of first (and second) generation fermions are very heavy [HI IS]- However, we 
do not assume a particular scenario for supersymmetry breaking (presumably described 
by high-scale dynamics), instead simply treating the relevant weak-scale soft breaking 
parameters as independent inputs. 

The MSSM contains four potential WIMPs, the three sneutrinos Vi and the lightest 
neutralino Xi- However, sneutrinos have full weak strength vector couplings to the Z, 
leading to a cross section for scattering on protons that is well above experimental con- 
straints, unless their masses are so large that their predicted relic density is well above 
the observed value. This leaves the lightest neutralino as unique WIMP candidate in the 
MSSM [ID]. 

While there have been many analyses of Xi as WIMP candidate, nearly all of them 
compute the relevant annihilation cross sections only at tree level. However, in order 
to achieve percent level precision, at least the leading radiative corrections need to be 
taken into account. These include QCD corrections for strongly interacting final states 
[TTJ H2J [E2J dH dSj, but at least leading electroweak corrections should also be included; 
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with "leading" we mean corrections that are parametrically larger than the naive one-loop 
factor a /it. 

Of course, a complete one-loop calculation would be ideal, and has in fact already 
been performed for a few points of parameter space [161 El HB]- However, calculating 
cross sections at one loop order involves computing a large number of diagrams. This not 
only needs a large effort even in the days of largely automatized calculation of one-loop 
amplitudes, the resulting expressions are also very complex, and their evaluation is very 
CPU intensive. With present-day computers the results of full one loop calculations can 
therefore not be used if the relic density has to be evaluated for many different choices 
of parameters. 

It is therefore useful to develop alternative ways to estimate the bulk of electroweak 
radiative corrections to the relic density. One class of diagrams that can lead to large 
corrections involves the exchange of a relatively light boson between the annihilating 
WIMPs. These corrections can even become non-perturbative if the ratio of the boson 
and WIMP masses is smaller than the relevant squared coupling a [191 EDI EH E2j , and 
they can remain significant even for larger mass of the exchanged boson. In the context of 
the MSSM these corrections are usually small [23], except in regions of parameter space 
where Xi is a heavy higgsino- or wino-like state fM\ I2"5] . 

In this article we are interested in another class of potentially large radiative correc- 
tions, which can be described by running or effective couplings. For example, it has been 
recognized some time ago [11] that running quark masses, or Yukawa couplings, evaluated 
at a scale near m^o, should be used in the calculation of the Xi renc density. Similarly, 
most calculations of the Xi renc density employ running electroweak couplings, based on 
the QED coupling at the weak scale (often, at scale Mz), rather than the on-shell value 
of a cm . 

In the limit of exact supersymmetry these electroweak gauge couplings not only ap- 
pear in Xi annihilation into final states containing gauge bosons, and in s— channel 
Z— exchange diagrams; they also determine the couplings of Xi to a Higgs boson and 
a neutralino or chargino, as well as the gauge contributions to the Xi couplings to a 
sfermion and a fermion, which have the structure of Yukawa couplings. However, when 
supersymmetry is broken, SM particles and their superpartners (for example, q and q) no 
longer have the same masses, which leads to different contributions to gauge and gaug- 
ino two-point functions. In particular, for energies between the squark and quark mass 
scales, the (s) quark superfields no longer contribute to the running of the gauge Yukawa 
couplings, but the quarks still contribute to the running of the true gauge couplings. 
Corrections to the gauge Yukawa couplings are the supersymmetric analog of the oblique 
corrections, and have thus been dubbed "super-oblique" corrections [261 127]P1 

While the first analyses of these super-oblique corrections only included terms that 
grow logarithmically with the masses of the heaviest superpartners [261 E2], references 
[301 l3Tj include the full quark and squark corrections to leptonic processes, thereby also 
including terms that approach a constant for large squark mass. Finally, in [32] this for- 
malism was extended to include all (s) fermion loop corrections involving gauge Yukawa 

§This apparent hard breaking of supersymmetry caused at one-loop order by soft breaking terms has 
first been noticed in the context of strong interactions in ref . |28j . 
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couplings to the chargino and neutralino two-point functions. This allowed the intro- 
duction of effective neutralino and chargino couplings to sfermions and fermions; here we 
apply the same formalism also to the gauge Yukawa couplings involving a Higgs boson 
and two neutralinos or charginos. 

Note that these corrections are enhanced not only by the logarithm of the ratio of the 
mass of the heaviest sparticle to the LSP mass, but also by the large number of fields 
contributing with the same sign (since the contributions to diagonal gaugino two-point 
functions go like the square of the relevant charge or coupling). The latter factor is more 
important, unless one is considering a scenario where some sfermions are more than an 
order of magnitude heavier than the LSP. A calculation [33] of radiative corrections to 
X2 Xi e+e ~ decays found that the bulk of the corrections to the partial width could 
be described by these effective couplings even in a scenario where all superparticles are 
relatively light. 

Since we include (certain) one-loop corrections to chargino and neutralino two-point 
functions, we also use one-loop renormalized neutralino and chargino masses. We use a 
recently developed [31] variant of the on-shell renormalization scheme [35] , where correc- 
tions to physical masses are small for (nearly) all combinations of parameters, in contrast 
to earlier schemes, which lead to large corrections when the SU (2) gaugino and higgsino 
mass parameters become close in magnitude. The effect of this mass renormalization on 
the relic density is therefore usually (although not always) small. 

This paper is organized as follows: in section [2] we outline our formalism and discuss 
the structure of the corrections included in the effective couplings. Section [3] contains 
technical details of our numerical calculation. Section H] is devoted to the discussion of 
numerical results and we conclude in section [5j Details of the renormalization scheme we 
chose are given in the Appendices. 



2 Formalism 



In this section we describe the formalism of effective couplings, closely following [32] , 
and discuss in some detail the nature of the corrections treated in this formalism. Here 
we are interested in corrections to the Yukawa couplings of the electroweak neutralinos 
and charginos. The required renormalization of the chargino, neutralino, Higgs and 
electroweak gauge sectors is described in the Appendices. 

At tree level the neutralino and chargino couplings in question can be written as 



N* al tan 9 W ( Q f - T 3 / ))P L f* L - tan 9 W Q f N al P R f* R f 



9 m f .~o 

Act 



c 



xt du 



V2M w f f ((3) 
gm u 



K h p L r R +N ahf p R r L f 



h.c. 



gmd 



V2M W sin 



V2M w cos(3 
xtV i2 P R udl + h.c 



xtU* 2 P L ud R 



(1) 



(2) 



The unitary 4x4 matrix N appearing in eq.Q diagonalizes the neutralino mass matrix, 
while the unitary 2x2 matrices U, V appearing in eq.([2]) diagonalize the chargino mass 



4 



matrix. Furthermore, g is the 577(2) gauge coupling, 9w is the weak mixing angle, T 3 
and Qf are the third component of the weak isospin and electric charge of fermion /, 
respectively, Pr^l — (1 ± 7s)/2 are the chiral projectors, m/ and are the masses of 
fermion / and of the W boson, respectively, and tan /3 is the ratio of vacuum expectation 
values (vevs) of the two neutral Higgs fields required in the MSSM . Finally, the function 
ff((3) appearing in eq.Q is sin f3 for up-type quarks, and cos (3 for down-type quarks and 
charged leptons; similarly, the index hf = 4 (3) for up-type quarks (down-type quarks 
and charged leptons). The chargino couplings to a neutrino and a charged slepton have 
the form of eq.(T5]) with u — > u, d — > £ (and m„ — > m v = in our treatment). The 
chargino couplings to a down-type quark and an up-type squark again have the form of 
eq.©, with u •<->■ d, U ■<->■ V and xt {xt) C \ where the superscript C denotes charge 
conjugation. 

Equations ([I]) and (j2J) have been written for neutralino and chargino mass eigenstates, 
but fermion and sfermion current eigenstates. In our numerical treatment we have ignored 
flavor mixing in both the fermion and sfermion sectors. Quark mixing is known to be quite 
small, and neutrino mixing is irrelevant for our purposes (since neutrinos remain very 
nearly degenerate and massless, compared to the sparticle masses). Moreover, sfermion 
flavor mixing is strongly constrained by experimental upper bounds on flavor changing 
neutral currents. However, in our numerical treatment we include fx, — fa mixing for 
third generation sfermions. 

Note that the first line in eq.([T]) and the first term in eq.([2]) come from gauge inter- 
actions, whereas the second line in eq.Q as well as the second and third term in eq.(j2D, 
which are proportional to SM fermion masses, originate from superpotential couplings. 
Technically the corrections we are interested in originate from fermion-sfermion loop 
corrections to the chargino and neutralino two-point functions, plus the counterterms 
from (s)fermion loops needed to render them finite. Note that the sum of these explicit 
two-point function corrections and counterterms has to be finite, since these are the 
only one-loop corrections to xff vertices involving (s)fermions /', /' of a different SU(2) 
multiplet. As usual, the finiteness of the results is an important check on the calculation. 

In [32] all these contributions were absorbed in corrections to the matrices N, U and 

V: 

U=U + AU, V = V + AV, N = N + AN. (3) 

Note that AU, AV, AN are not counter terms; rather, they describe finite corrections 
to the couplings of charginos and neutralinos to sfermions and fermions, and to Higgs 
bosonsEl In our renormalization scheme, as in ref.[32], the mixing matrices are not renor- 
malized at all, since we allow off-diagonal wave function renormalization counter terms. 

*A perhaps clearer, but bulkier, notation would be A(gU), A(gV), A(gN), which can be obtained 
from eqs.@ by multiplying both sides with g; note that even the contributions to the chargino and 
neutralino couplings ((lj and (0) that originate from superpotential couplings have a factor of g in front. 
For easier comparison we here stick to the notation introduced by Guasch et al. [52] . 
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Explicitly, the corrections defining the effective couplings are given by 




+ V 3 2^Z± 



■H 

Pa J 



(4) 



In the first four equations, for AU and AV, the subscript j in the last term must always be 
different from i appearing on the left-hand sides, i.e. these are contributions from purely 
off-diagonal two-point functions, just like the corresponding contributions to the AiV. 
The counterterms Sg and S tan 8w are given in Appendix A, while 5 sin (3 and S cos are 
given in Appendix C. Finally, the diagonal and off-diagonal wave function counterterms, 
5Zi and 5Zij, are given at the end of Appendix D; here the superscript L and R refer to 
chirality, while the superscript + refers to charginos. 

We emphasize that AN, AU and AV, as described above, are process independent in 
the sense that they do not depend on the details of the external state [the flavor of the 
(s)fermion or the identity of the Higgs boson]. 

An important feature of these corrections is that heavy sparticles do not decouple; to 
the contrary, the size of the corrections increases logarithmically if some sfermions are 
much heavier than the charginos or neutralinos whose couplings are being computed. This 
observation is especially relevant in the light of recent LHC data which imply stringent 
lower bounds on the masses of first and second generation squarks, but allow relatively 
light superparticles with only electroweak interactions. As an illustration let us first 
review the origin of these non-decoupling effects only including gauge interactions. 

Unbroken supersymmetry implies that the fermion-fermion-gauge boson (gauge) cou- 
pling g is the same as the corresponding sfermion-fermion-gaugino coupling g. Since 
supersymmetry is broken at a high scale (compared to the light fermion masses), the 
heavier sfermions and fermions will decouple from the running of g below the sfermion 
mass scale m^. However, the corresponding fermion will still contribute to the running 
of g. Therefore, g and g "run apart" below the highest supersymmetry breaking scale. 



6 



We use the running coupling a(Mz) as tree-level input coupling, which is defined at 
the fixed scale Mz- Consider a spectrum where sleptons, charginos and neutralinos have 
similar masses to each other, while squarks have much higher masses, mj ~ m^. <C rriq. 
We are interested in processes at the slepton or neutralino mass scale. On the other hand, 
full supersymmetry is attained only at scale rriq. Since between scales rriy, and rriq the 
leptonic sector is still supersymmetric, leptons and sleptons contribute in the same way 
to the running of g and g in that range. Altogether we thus have: 

9{ m x) = 9i m q) -Pejlog — 

m x 

= g(mq) - fi t flog— 1 
m x 

= g{M z ) + ft, log ^ + /^log ^ - ^-log ^ 

77? - 777 - 

= 0(M Z ) + P e , q log + log -1. (5) 

Mz rn^ 

Here the particles contributing to the gauge beta function have been listed explicity as 
subscripts; in the last step we have used the fact that to one-loop order the contributions 
of different (s)particles simply add up, i.e. (3 £q i = (3 t i+/3 q . Since SM fermions contribute 
with positive sign to the gauge beta functions, both logarithmic factors in the last line of 
eq.© are positive [26| 127]. 

We use on-shell renormalization both in the electroweak and in the x sector. Our 
effective x% couplings can thus be understood as running couplings at scale m^. This 
is the appropriate energy scale to describe x% decays, where the effective couplings have 
been shown to reproduce the bulk of electroweak radiative corrections [33], as well as for 
the annihilation of two neutralinos, which is the main topic of the present investigation. 
It may not be the best choice for sfermion decays, where these couplings have first been 
introduced [32], nor for slepton production at a lepton collider via x exchange, where the 
X is always quite far off-shell. 

The existence of Higgs superfields, whose superpotential couplings give rise to matter 
fermion masses and whose fermionic members mix with the gauginos, considerably com- 
plicates the picture. Fortunately the chirality structure allows to distinguish "gauge" and 
"superpotential" contributions to the xff couplings: the former couple fz to fz and - 
for neutralinos coupling via hypercharge - fjt to /«, whereas the latter couple fz to fa 
and f R to fz- 

At one-loop level true Yukawa couplings A do not renormalize true gauge couplings 
g. The above discussion of the 0(g 3 ) corrections to the "gauge" contribution of xff 
couplings therefore indicates that there are no 0(g\ 2 ) corrections to these couplings. This 
is in fact true in the absence of gaugino-higgsino mixing, or, to good approximation, for 
gaugino-like states, where 0(g\ 2 ) corrections are suppressed by two factors of higgsino- 
gaugino mixing angles. 

However, the "gauge" couplings of higgsino-like states do receive corrections 0(g\ 2 ), 
where (s)top, (s)bottom and (s)tau loops all can contribute. These corrections are sup- 
pressed by one factor of a higgsino-gaugino mixing angle e; however, since this is true 
already for this tree-level coupling, the relative size of the correction is not suppressed. 
Technically, these corrections involve the "vector" parts of diagonal as well as off-diagonal 
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X two point functions, as well as the "scalar" parts of off-diagonal two-point functions; 
see eq. f[26l) for the decomposition of fermionic two-point functions. In the latter case, the 
product gXfirif has been counted as O(Xje), where the factor e describing suppression by 
higgsino-gaugino mixing arises from the combination of the factor My/ in mj oc XfMw/g 
and the factors in the expressions for the off-diagonal wave function counterterms 
given in eqs. (l38i) and (1391) . 

Since there are no unsuppressed 0(gX 2 ) corrections to the "gauge" couplings of 
gaugino-like states, these corrections cannot be understood in terms of running g cou- 
plings. Instead, in the framework of the effective theory defined by integrating out heavy 
sfermions, they can be interpreted as contributions to the running of those (off-diagonal) 
elements in the chargino and neutralino mass matrices that mix higgsino- and gaugino- 
like states. 

This interpretation is supported by the observation that the "superpotential" cou- 
plings A of gaugino-like states receive one-loop corrections 0(Xg 2 ), 0(X 3 ), where again 
all heavy (s)fermions contribute, that are only linearly suppressed by gaugino-higgsino 
mixing. Since these tree-level couplings suffer a similar suppression, the relative size of 
the corrections is again not suppressed. In contrast, the "superpotential" coupling of 
higgsino-like states only receives unsuppressed 0(X 3 ) corrections, while 0(Xg 2 ) correc- 
tions are suppressed by two powers of higgsino-gaugino mixing angles; moreover, the for 
our application most relevant coupling oc A r in this case does not receive unsuppressed 
corrections from the (s)top sector. This is consistent with the observation that only true 
Yukawa couplings contribute to the beta functions of other Yukawa couplings, if one 
restricts oneself to matter (s)fermion loops; since matter loops contribute with positive 
sign to the beta functions of superpotential couplings, the corresponding corrections are 
positive, just as in the pure gauge sector 

Finally, the "gauge" couplings of higgsino-like states also receive 0(g 3 ) corrections, 
again suppressed by one factor of higgsino-gaugino mixing; these originate both from 
the running of the "gauge" coupling g and from 0(g 2 ) corrections to gaugino-higgsino 
mixing. Similarly, the "superpotential" couplings A of gaugino-like states receive 0(X 3 ) 
corrections, suppressed by one factor of higgsino-gaugino mixing, that can be understood 
from the running of A. 

Altogether, the tree-level and logarithmically enhanced yii couplings, where I stands 
for a lepton, can be summarized by table [TJ where we have not distinguished between 
g and g, nor between the superpotential coupling A and the corresponding contribution 
A to xff couplings. Of course, our effective couplings also include terms that are not 
enhanced logarithmically for large sfermion mass. The corrections that vanish as some 
inverse power of the mass of the sfermion in the loop do not follow the pattern of this table; 
for example, there are 0(gX 2 m 2 ^/m 2 ) corrections to the "gauge" couplings of gaugino- 
like x states. However, these decoupling corrections are usually small in our numerical 
applications. 

We repeat that the non-decoupling corrections from the gauge sector are enhanced 
by multiplicity factors in addition to possibly large logarithmic factors. Moreover, the 

^In the full theory there are 0(g 2 X) corrections to the beta function of A, but these are due to diagrams 
where members of gauge superfields run in the loop. These contributions cannot be included consistently 
into our effective couplings. 
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coupling 


tree-level 


X — higgsino 

one-loop 


X - 
tree-level 


n gaugino 

one-loop 


XhtR, X^rIl 


0{eg) 
0(Xe) 


<D(eg\ egX 2 ) 
0(\l XAl e 2 A^ 2 , e 2 A £ A?) 


0(g) 
0(eX e ) 


0{g\ e 2 gX 2 ) 
0(eX e g 2 , eX e X 2 ) 



Table 1: Comparison of tree-level and one-loop contributions to slepton-lepton currents 
of different chiral structure, for higgsino- and gaugino-like charginos and neutralinos X- 
Here g stands for an electroweak gauge coupling, Xi for the Yukawa coupling of lepton I, 
and e for one factor of gaugino-higgsino mixing, which is assumed to be small here. A 
coupling A without subscript means that all superpotential couplings contribute without 
additional suppression. We do not distinguish between g and g, nor between A and A, in 
this table. 

effective couplings capture all corrections oc X T Xf, X r Xf to tree-level couplings oc A r . They 
do, however, not capture all corrections oc A^ to tree-level couplings oc X t ; since the latter 
also receive QCD corrections, which are also not included in our effective couplings, their 
application to xqq couplings involving superpotential couplings is probably of limited 
usefulness. 

As mentioned above, the corrections described by eqs.fll]) can all be computed from 
two-point functions; see the Appendices for details. At a similar level of effort one can 
include corrections to the (pole) masses of the charginos and neutralinos. Of course, these 
are always computed from two-point functions only; we can, and do, therefore include 
the complete correction to the masses, including corrections from the gauge and Higgs 
sectors in addition to those from the matter (s)fermion sectorQ We do so using a recently 
introduced [31] variant of the on-shell scheme, where we consider the masses of the (more) 
wino-like chargino and of the the (most) bino- and higgsino-like neutralinos as inputs, 
which are not altered by the corrections. In ref. [51] it was shown that the counterterms 
to M\, M 2 , and \x then remain well behaved in most of the parameter space, in contrast 
other versions of the on-shell scheme [35]. As a result AN, AU and AV also remain well 
behaved even when M 2 = 

The expressions 01]) were originally introduced [32] as corrections to couplings of 
charginos and neutralinos to fermions and sfermions. However, in addition, the correc- 
tions AN a i, AN a2 , AUn and AVu can be used to (partly) absorb radiative corrections 
to all the vertices which involve a gaugino component of a chargino or neutralino but 
do not include a gauge boson. These include couplings of neutralinos and charginos to 
Higgs bosons, the gaugino component of which we also modify. However, the corrections 
AiV Q 3, AiV Q 4, AUi2, AVi2 should not be applied to these Higgs vertices, but only to ver- 
tices involving a (massive) fermion and one of its superpartners; as discussed above, part 
of the latter corrections can be understood as process independent corrections to chargino 
and neutralino couplings originating from the superpotential, whereas the couplings of 

* Computing the corrections to neutralino and chargino masses entails again computation of the coun- 
terterms 5M\y, Sg, 6cos9w and <5cos/3, now including all one-loop corrections. In the calculations of 
the one-loop corrected masses these counterterms will therefore have different numerical values than 
those appearing in eqs.@, where only corrections from (s)fermion loops must be included in order to 
ensure finiteness of the effective couplings. 
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Higgs bosons to neutralinos and charginos are pure gauge couplings in the limit of exact 
supersymmetryJl] Similarly, a vertex with a gauge boson involves a true gauge coupling. 
The running of the corresponding gauge coupling (the U(l)y and SUl(2) gauge couplings 
in our context) takes care of all process independent radiative corrections to these ver- 
tices, which should therefore be written in terms of the original U, V, N matrices, without 
the corrections of eqs.(|4]). 

Before presenting numerical results for corrections to neutralino annihilation cross 
sections, we wish to show some results for the corrections to neutralino couplings, i.e. the 
AN a /3. In Fig. [H we plot the relative size of these corrections, AN a p/N a p, for selected 
values of a and /3, against a common (i.e. generation and flavor independent) soft SUSY 
breaking sfermion mass Msusy- We assume the following parameters (all mass parameters 
are in GeV): 

tan/3 = 10, M A = 500, M x = 100, M 2 = 300, M 3 = 1200, \i = 600, A = 0. (6) 

Since fi — M\^, Mz, gaugino-higgsino mixing is suppressed, so the results of table [1] 
are applicable. Specifically, Xi is bino-like, x° an d Xi are wino-like, and xj], Xa an d xt 
are higgsino-like. 




Figure 1: SN a/3 = , as given by eq.(jl]) against Msusy, where MgusY denotes the 

N a /3 

common soft SUSY breaking sfermion mass. 

The left frame shows corrections to couplings that are not suppressed by gaugino- 
higgsino mixing. Recall that the subscript a in N a p refers to the mass eigenstate (ordered 
by increasing mass), while j3 = 1, 2, 3, 4 refers to the bino, wino, down-higgsino and up- 
higgsino component, respectively. Clearly the biggest correction are those to N44, which 
involve superpotential couplings of the "up-type" higgsino, in particular the top coupling 
X t . There are unsuppressed C(A^) contributions to ^iV^; this case is analogous to the 
0(\f) corrections to the x£l^r couplings in table [lj For the given, rather small, value 

t Another way to see this is to note that it makes little sense to include higgsino two-point function 
corrections to a Higgs-higgsino-gaugino vertex, without including closely related two-point function 
corrections on the Higgs boson line. 



10 



of tan/3, 5N 33 is far smaller, since it only receives contributions oc Xf, X 2 T) plus C(e 2 ) 
contributions. Not surprisingly, if tan/3 ~ mt/m&, SN33 is of similar size as &/V44. 

Among the corrections to the dominant couplings of the gaugino-like neutralinos, SN22 
increases significantly faster with increasing Msusy than 8Nn does simply because the 
SU{2) gauge coupling is larger than the U(l)y coupling; this is only partly compensated 
by the fact that SU{2) singlet (s)fermions do contribute to 5Nu. 

The right frame of Fig. [TJ shows corrections to couplings that are 0(e) at tree-level. 
The biggest corrections are to N u and N 2 4, involving "superpotential" couplings of 
gaugino-like states. These corrections are dominated by O(Xf) running coupling ef- 
fects, and are in fact of similar magnitude as shown in the left frame; recall that 
we are showing the relative size of the one-loop correction here, so that the 0(e) factor 
appearing in all couplings shown in the right frame cancel out. 

The corrections to the "gauge" couplings of the higgsino-like states and are 
much smaller. These receive running coupling contributions proportional to the squared 
U(l)y and SU{2) gauge coupling, respectively, but are much smaller than 5Nn and 5N 2 2 
shown in the left frame. This indicates that the corrections to gaugino-higgsino mixing, 
which are dominated by O(Xf) contributions, tend to reduce the couplings in this case. 
In fact, the total coefficient of the correction oc log(MsusY) to 5N^i is negative. The 
fact that the overall correction nevertheless remains positive even at Msusy = 30 TeV 
illustrates that the terms that approach constants for large sfermion masses can also be 
very important. 

This completes our discussion of the effective couplings. Before presenting numerical 
results for the annihilation cross section and relic density of the lightest neutralino, we 
briefly describe the implementation of the effective couplings using publicly available 
software. 

3 Implementation 

In order to predict the relic density of the lightest neutralino at one-loop level, we have 
modified one of the publicly available codes for relic density calculation, micrOMEGAs 
(version 2.2.CPC.i) [36]. The algorithm describing the flow of relevant parameters is 
sketched in Fig. |2j 

We use SuSpect [37] to compute the "tree-level" spectrum of superparticles and 
Higgs bosons. SuSpect actually includes various loop corrections in converting input 
DR parameters into on-shell parameters and, hence, physical masses. However, in the 
neutralino-chargino sector these corrections are complete only in the absence of higgsino- 
gaugino mixing. Since we implement full one-loop corrections to all x masses in the 
on-shell scheme, we had to slightly modify [31] the x mass matrices used in SuSpect, 
using on-shell W and Z masses and the on-shell definition of the weak mixing angle 6w 

The output generated by SuSpect (in the format of the Supersymmetric Les Houches 
Accord, SLHA [SB]) is used as input of micrOMEGAs. This program first computes two 
external libraries, written by us, which calculate the effective couplings and the renormal- 
ized masses, respectively. We have used FeynArts, FormCalc [391I10JUT] and LoopTools 
[12] to compute the relevant diagrams and counterterms. Feynman gauge has been used 
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micrOMEGAs 
(calls the renormalized mass 
library and the effective 
couplings library. Redefines the 
chargino and neutralino masses.) 
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Effective couplings 









CalcHEP 
(with modified vertices wherever 
applicable, reads AN, AU 
and AV from a file generated 
by effective couplings library) 



relic density with 
desired radiative corrections 



Figure 2: Flowchart describing the interaction of the publicly available codes SuSpect and 
micrOMEGAs with our subroutines calculating corrections to x masses and couplings. 



in the calculation of one-loop corrections to \ masses; the effective couplings do not 
include loops involving gauge bosons. For regularization we have used the constrained 
differential renormalization method jl3]. At one-loop level this method has been proven 
to be equivalent [H] to regularization by the dimensional reduction method [H]. Note 
that, unlike the effective couplings, the renormalized masses involve contribution from all 
MSSM (s)particles; the two external libraries therefore use independent calculations of 
the relevant counterterms, as further described in the Appendices. In the loop diagrams, 
we for simplicity used running DR sfermion masses defined at the electroweak scale, as 
given by SuSpect the use of a different renormalization scheme for the sfermion masses 
inside the loops affects our results only at two-loop level. Moreover, we ignore L — R 
mixing for first and second generation sfermions, and do not include any flavor mixing. 

MicrOMEGAs calculates the Xi (co-) annihilation cross sections with the help of CalcHEP 
[15] . We modified CalcHEP to read the effective couplings AU, AV and AN. These 
are used in all xf f couplings; the corrections to gaugino couplings AUn, AVu, AN a i 
and AN a2 are also used in Hxx vertices, but the corrections to higgsino couplings 
AU a , AVi2, AN a3 and AN ai are not used in these vertices, and the couplings between x 
states and gauge bosons are not modified at all. 

We already mentioned that we use on-shell renormalization in the electroweak sec- 
tor. This refers to the W and Z masses, the weak mixing angle 8w as well as the 
running electromagnetic coupling a(Mz)- For consistency we thus have to use on-shell 

*The reason is that the version of SuSpect we were using implemented certain O(X^) corrections 
to stop masses, but not to sbottom masses, thereby creating spurious hard breaking of SU(2) in the 
stop-sbottom sector, with sizable effects on certain electroweak counterterms. 
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values for these quantities also in the tree-level calculation. This required additional 
changes of micrOMEGAs, which uses an on-shell Mz, but the MS value of a(Mz), which 
is nearly 1% larger than the on-shell value, and an effective weak mixing angle such that 
sin 2 9w ~ 0.231 (compared to 0.222 in the on-shell scheme). We use this modified version 
of micrOMEGAs only for the one-loop corrected cross section. Since micrOMEGAs currently 
is the "industry standard" for calculations of the x? re hc density, we use its original 
version when evaluating the final corrections to the (co-) annihilation cross sections and 
relic density. The relative corrections to the annihilation cross sections are thus defined 
as 

tree, on— shell _i_ A ^.1— loop -.tree, orig 

5a = , (7) 

Q-tvee, orig ' v ' 

where <7 tree ' ong is the prediction of the original version of micrOMEGAs using MS values of 
a (M z ) and 6 w \j\ We emphasize that using a tree ' on , computed with on-shell values for 
8\y and a(Mz), as reference value would have led to significantly larger loop corrections 
8a for bino-like xli largely because a(Mz) is smaller in the on-shell scheme, as mentioned 
above. Note finally that we use on-shell values of aiMz) and 9w only in those vertices 
which we correct by the effective couplings. In all other vertices we continue to use 
micrOMEGAs default values of parameters; this includes all couplings of gauge bosons, 
as well as couplings of Higgs bosons to SM particles. This ensures that 5a defined in 
eq.© only shows effects of loop corrections we include but standard micrOMEGAs doesn't 
included 



4 Numerical Results 

Before presenting numerical results for a benchmark scenario, a qualitative discussion of 
the usefulness of the effective couplings for the calculation of the Xi re bc density might 
be in order. 

We expect the effective couplings to be most useful for dominantly bino-like neutralino 
|18j . A bino-like neutralino dominantly annihilates through t and u channel sfermion 
exchange to a pair of (lighter) fermions. Effective couplings absorb significant radiative 
corrections in both vertices. These will likely dominate the total electroweak corrections 
whenever squarks are significantly heavier than neutralinos and sleptons, since this leads 
to logarithmic enhancement of our corrections, as discussed above. The most important 
final states will then be £ + £~ pairs. Since the cross section is oc |iVii| 4 , the relative size 
of the correction can be estimated as 45 Nn. 

However, given existing lower bounds on slepton masses [IS], pure slepton exchange 
usually gives too small an annihilation cross section, i.e. too large a Xi re li c density (in 
standard cosmology). This can be corrected through co-annihilation with a slepton i, 
which also involves the Xi^ vertex. 

^Recall, however, that we slightly changed the off-diagonal entries in the chargino and neutralino 
mass matrices in SuSpect; the resulting x masses and mixing angles are also used by the "original" 
micrOMEGAs in our calculation. 

••■For example, changing a(Mz) and 6w everywhere would also lead to changes of the Higgs spectrum, 
which can lead to large changes of the annihilation cross section near Higgs poles; these would be entirely 
due to changes of input parameters, not due to radiative corrections. 
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The annihilation cross section of bino-like neutralinos can also be enhanced sufficiently 
(or even by too much) if there is a neutral Higgs boson with mass near twice the neutralino 
mass, the CP-odd state A being most efficient for this purpose. We also use an effective 
coupling for the gaugino components in the XX Higgs vertices. The flavor independent 
contributions to the two point function corrections to the gaugino components remain 
finite, and can absorb significant corrections to these processes. 

Given lower bounds on chargino masses, higgsino- or wino-like neutralinos will an- 
nihilate dominantly into gauge boson pairs W + W~, ZZ. These cross sections involve 
true gauge interactions, i.e. logarithmically enhanced loop corrections can be absorbed 
into the running of the SU(2) [and, in case of higgsinos, U{l)y] gauge coupling. Our 
effective couplings can still capture leading corrections to subleading // final states for 
the annihilation of wino-like neutralinos; they are also likely to be quite important in the 
vicinity of possible s— channel Higgs resonances. 

Let us now turn to a discussion of numerical results, based on the following benchmark 
scenario: 

rriq = 1.5 TeV, M 3 = 1.2 TeV, m~ lh = 0.55 TeV, m~ lR = 0.50 TeV, 
H = 0.6 TeV, m A = 0.5 TeV, A 3 = 1.0 TeV, tan(/3)(M z ) = 10 . (8) 

All parameters except the on-shell mass of the CP-odd Higgs boson are running 
DR parameters, in most cases defined at the scale of electroweak symmetry breaking 
(EWSB), which in SuSpect is defined as y/fn^ m~ t2 ~ 1.5 TeV in our benchmark scenario. 
This choice minimizes one-loop corrections to the effective potential for the neutral Higgs 
fields jlTJ 08]. Different values will be chosen for the SU(2) and U(1)y gaugino masses 
M 2 and Mi in subsequent figures, as stated below. 

The input value of tan/3 is defined at scale M z , which is the SuSpect default. We 
have used tan/3 DR at the EWSB scale, approximately equal to 9.6, for our calculations. 
We do not include any plots against tan /3, since neutralino masses and mixing matrices 
are not very sensitive to it. 

In all parameter scans we ensure that x\ is the lightest supersymmetric particle. 
Since we assume exact R parity, \\ remains a dark matter candidate in these regions of 
parameter space. However, since the main goal of this paper is to show the relative size of 
the radiative corrections captured by effective couplings, we do not restrict ourselves to 
regions of parameter space giving the correct relic density in standard cosmology. Finally, 
in all figures showing x\ annihilation cross-sections, these have been calculated assuming 
the relative velocity between the annihilating neutralinos in their center of mass frame 
to be v = 0.3, which is typical for the epoch when the Xi decoupled from SM particles. 

4.1 Mi dependence 

In Fig. |3t we plot ANu, the correction to the effective bino coupling of Xi, against Mi, for 
M 2 = 400 GeV. Note that we show the total correction here, not the relative correction 
shown in Fig. [TJ For Mi < 410 GeV, Xi is bino-like and ANu — 0.018 is essentially 
independent of M\. In this region of parameter space it is dominated by RG-like 0(g Y ) 
contributions, where gy is the U(i) Y gauge coupling, and is thus positive. For M 1 > 410 
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Figure 3: Effect of radiative corrections as function of the bino mass parameter Mi (in 
the DR scheme at the electroweak symmetry breaking scale), for our benchmark point 
(JSJ) with M 2 = 0.4 TeV. We show the effective coupling Nu defined in eqs.fllj) (top left); 
the relative correction to the predicted Xi relic density (top right); the relative correction 
to (t(xiXi ~^ e + e~) (bottom left); and the relative correction to the total x\ annihilation 
cross section (bottom right). The last two plots are for fixed relative velocity v = 0.3. 
All relative corrections have been computed as in eq.flZJ). 



GeV, Xi becomes wino-like. This has little influence on the relative size of the correction 
to the bino coupling, but, since |iVn| becomes much smaller than in the region of bino- 
like Xv AiVn is also greatly reduced. Correspondingly, in this part of parameter space 
AiVn no longer gives the most important correction to Xi annihilation. 

The relative size of the correction to the xi P a h annihilation cross section into an e + e~ 
pair, defined as in eq.(JZJ), is shown in Fig. [3b. Over most of parameter space this cross 
section receives its dominant contributions from selectron (in particular, ejt, due to the 
larger hypercharge) exchange in the t— and u— channels. For Mi < 400 GeV, where Xi is 
bino-like, the correction is then essentially given by 4AAT n ; the relative correction shown 
is smaller than this, because our reference cross section does not use on-shell values for 
a(Mz) and $w, as explained in eq.(|ZI). This indicates that it is indeed advantageous to 
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use the MS value of aiMz) if one is performing a tree-level calculation, since it reduces 
the size of the radiative corrections^ 

An exception occurs for M\ ~ 47 GeV, where Xi annihilation can occur through the 
exchange of an almost on-shell Z boson. The very large enhancement factor (Mz/Tz) 2 — 
1300 overcompensates the coupling suppression |A^ 3 | 2 — | A^^j 2 of this contribution to the 
matrix element. Since we do not modify the ZxiXi coupling, the size of the correction 
decreases at the Z pole. 

As noted above, for M\ > 410 GeV Xi becomes wino-like. Here the corrections, which 
are now dominated by terms oc ANu, become much larger; note that the cross section 
is now dominated by e~i exchange. We will see in Fig. |4^l that | <JiVi2 1 — 1.3% increases 
the absolute size of the wino coupling, again due to RG-like effects. Moreover, recall 
that our tree-level calculation uses on-shell values of both aiMz) and 9w As a result, 
our tree-level SU{2) gauge coupling is about 1.5% larger than that used by micrOMEGAs. 
Altogether our e.L&x\ wino coupling is therefore about 3% larger than in micrOMEGAs, 
leading to an about 12% enhancement of the Xi annihilation cross section into e + e~ 
pairs in the deep wino-like region. Note that for wino-like Xi> our choice of tree-level 
couplings leads to smaller radiative corrections than the choice made in micrOMEGAs. 
This is in contrast to experience with LEP observables, where the "effective" or 
MS [50] 6vM leads to smaller radiative corrections than the electroweak on-shell scheme. 
Recall, however, that for bino-like Xi the use of MS values for 9w and a(Mz) reduces 
the size of radiative corrections. 

The relative correction to the total Xi annihilation cross section, summed over all 
kinematically accessible final states, is shown in Fig. [3H. It exhibits a considerably more 
complicated dependence on M\ than the cross section for the e + e~ final state does. The 
reason is that different final states become important for different ranges of M\. 

Near the Z pole the total annihilation cross section is again dominated by Z exchange, 
and the correction is small. For slightly larger Mi, the s— channel exchange of the lighter 
CP-even Higgs state h becomes (nearly) resonant. As noted earlier, we correct the 
gaugino component of the hXiXi coupling. Since the h exchange matrix element is linear 
in this coupling, the correction is considerably smaller here than in regions of parameter 
space where the total annihilation cross section is dominated by I exchange in the t— 
or u — channel. The correction therefore increases in size beyond the h— pole, reaching a 
maximum close to the relative size of the correction to the cross section for the e + e~ final 
state at M x ~ 80 GeV. 

At larger values of M\ additional channels begin to become accessible: annihilation 
into W + W~, ZZ, Zh and hh final states becomes possible at M 1 ~ 84 GeV, 94 GeV, 
106 GeV and 118 GeV, respectively. Since the gauge vertices are not corrected and 
in many diagrams the Higgs vertex again only occurs linearly in the matrix element, 
the corrections begin to decrease; however, the total annihilation cross section is still 
dominated by final states. 

*The correction is further reduced because the small wino coupling of Xxi whose magnitude is also 
increased by the correction AiVfo, reduces the overall e^ex^ coupling. However, this effect is numerically 
not very important, since the wino coupling remains very small, and the exchange contribution is 
about 20 times smaller than the in exchange contribution; note that these contributions do not interfere 
in the limit m e — > 0. 

^These two definitions of the weak mixing angle differ conceptually, but are numerically very similar. 
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At M\ ~ 176 GeV annihilation into it opens up. Since this cross section is not 
P— wave suppressed, it contributes significantly to the total annihilation cross section, 
although the relevant matrix element is suppressed either by the large stop masses (for 
the t— and u— channel diagrams) or the small X1X1 Higgs couplings (for the s— channel 
Higgs exchange diagrams). The opening of this channel leads to a further reduction of 
the relative size of the corrections. 

For yet larger values of Mi, s— channel exchange of the heavy Higgs bosons H and 
A begins to dominate the total cross section. Since the corresponding matrix elements 
are again linear in the coupling we modify, the relative size of the corrections continues 
to decrease until, near Mi = 300 GeV, in rapid succession the H ± W T , ZA, ZH and 
finally hA and hH final states become accessible. Many of these final states contribute 
significantly to the total annihilation cross section; note that, unlike hh and £ + £~ produc- 
tion, hA production is possible from an S— wave initial state, i.e. the cross section is not 
suppressed by a factor v 2 = 0.09 in our calculation. Moreover, final states with two Higgs 
bosons can be produced via the exchange of higgsino-like states in the t— or w— channel; 
the corresponding diagrams are not suppressed by small couplings, and the suppression 
0{M\lii) by the heavy higgsino propagator is not very severe for Mi ~ 300 GeV. Since 
these diagrams involve two XX Higgs couplings, the relative size of the corrections begins 
to increase again for Mi > 300 GeV. 

Finally, for Mi > 410 GeV the wino component of x? begins to grow quickly. This 
leads to a rapid increase of the annihilation cross section into W + W~ pairs, mostly via 
t— and u — channel Xi exchange diagrams. Since we do not modify any of the relevant 
couplings, the corrections quickly decrease, and remain very small once |i\T 12 | — 1. 

The relative size of the correction to the predicted Xi re li c density is shown in Fig. [3b. 
For Mi < 350 GeV it is essentially the mirror image of the correction to the total anni- 
hilation cross section, which we just discussed, because to good approximation the relic 
density is inversely proportional to the total annihilation cross section. There are some 
minor differences, because the evaluation of the relic density involves thermal averaging 
over a range of velocities, followed by integration over a range of temperatures, whereas 
the results of the cross section had been computed for a fixed relative Xi velocity, as men- 
tioned above. As a result, the relative corrections to the relic density are slightly smaller 
in magnitude than those to the total annihilation cross section. This can be understood 
from the usual (semi-) analytical solution [51] of the Boltzmann equation determining the 
relic density. For a P— wave initial state, the relic density is inversely proportional to the 
product of the square of the decoupling temperature and the annihilation cross section. 
Our positive corrections increase the annihilation cross section, but also slightly decrease 
the decoupling temperature; because the latter depends only logarithmically on the an- 
nihilation cross section, the predicted relic density is indeed reduced, but by a slightly 
smaller amount than if the decoupling temperature was held fixed. The rapid change of 
the total annihilation cross section, and resulting change of decoupling temperature, also 
leads to a small kink at M\ = 256 GeV in the relative corrections to the relic density, 
right at the H, A exchange poles. 

For Mi > 350 GeV, co-annihilation with the wino-like states x!> an d Xi begins to 
be important. For Mi > 410 GeV, Xi itself is wino-like. A small discontinuity occurs at 
the cross-over value of Mi, since in our on-shell renormalization of the x sector the mass 
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of the most bino-like state is an input, i.e. is not corrected, while the mass of the most 
wino-like neutralino does receive (small) corrections. For Mi > 410 GeV we therefore 
not only correct certain Xx couplings, but also its massU 

In this region, co-annihilation with the (now bino-like) X2 no longer plays much of a 
role, but co-annihilation with Xx remains very important. This can lead to W ± Z final 
states, which are not affected by our corrections, as well as //' final states. The latter can 
be produced both via exchange in the s— channel and by /, /' exchange in the t— or 
u— channel. Squark exchange diagrams are strongly suppressed by the large squark masses 
in our benchmark point. For leptonic final states we find strong cancellation between W 
and I exchange diagrams, which enhance the importance of the corrections to the Xx^ 
as well as xt^' vertices. However, since most relevant (co-) annihilation cross sections 
are basically not sensitive to our corrections, the total size of the corrections to the Xx 
relic density remains below 1% here. Note that the correction to the relic density here 
has the same sign as the correction to the Xi annihilation cross section shown in Fig. [3U, 
since co-annihilation processes are more important than annihilation reactions. 

Our corrections lead to an even more perfect cancellation between W ± and £ exchange 
contributions to the production of leptonic final states in x\Xx co-annihilation, thereby 
slightly increasing the predicted Xx Ie ^ c density. 

Recall that all quantities in Fig. [3] are plotted against the DR parameter M\. The 
conversion of this, and other input DR parameters in the chargino and neutralino mass 
matrices, into on-shell parameters made by SuSpect is exact only in the limit of vanishing 
mixing between the x states. According to ref.[52] this reproduces the exact on-shell 
masses quite accurately. This is consistent with our finding that one-loop corrections 
to x masses are small, at least in our version of the on-shell scheme. However, in the 
region of strong bino-wino mixing, where bino-wino co-annihilation is important, a small 
change of parameters in the x mass matrices can lead to relatively large changes of mixing 
angles. In particular, using running DR values of My/, Mz and 6y/ m the off-diagonal 
elements of the x mass matrices, as in the original SuSpect, leads to results for A^ 12 and 
N 2 2 which differ by up to 5% from those we obtain by using on-shell values of My/, Mz 
and By/ everywhere, whereas the values of the x masses agree to within a fraction of a 
percent between the two calculations. 

It might seem peculiar that the original SuSpect uses running, scale-dependent pa- 
rameters in what are meant to be on-shell mass matrices. However, while on-shell masses 
are well defined as the poles of the (real parts of the) corresponding propagators, there 
is no correspondingly simple definition of a (non-diagonal) "on-shell mass matrvtf . In- 
deed, we could have performed an on-shell renormalization of the x sector where the 
off-diagonal elements of the tree-level mass matrices were defined in terms of running 
My/, Mz and 9y/', of course, (the finite parts of) certain counterterms would then dif- 
fer from their values in our calculation. If one insists on using DR input parameters, 
as in the supersymmetric Les Houches accord the most consistent treatment is to 
use DR renormalization throughout. On the other hand, if one insists on using on-shell 
renormalization, as we (and most other calculations of quantum corrections in supersym- 
metric theories) do, it seems more consistent to use the masses of the three input states 
as free parameters, rather than some DR masses. If we compare results of our modified 

•'•The discontinuity in ANn shown in Fig. [3^ has the same origin. 
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SuSpect with the original one for fixed values of the three input masses (e.g. of xf, Xi 
and X3 f° r Mi < 410 GeV), rather than for fixed Mi, M2 and /1, the two calculations 
also give very similar x gaugino components of the two lightest neutralinosJ§ This does 
not solve the problem of fixing the definition of on-shell mass matrices starting from DR 
input parameters; this obverservation is nevertheless reassuring, since it indicates that 
the precise definition of the on-shell mass matrix may not matter very much once actual 
measurements are used to fix the parameters of the theory. 

Finally, we mention that the Xi reuc density attains the desired value Q^oh 2 ~ 0.11 
for Mi ~ 230 GeV, ~ 265 GeV (to either sides of the H, A poles) and ~ 375 GeV (where 
the wino component of Xi as wen as co-annihilation with the wino-like states x 2 an d 
Xf become significant). For 230 GeV < Mi < 265 GeV as well as for Mi > 380 GeV 
the LSP could contribute a sub-dominant component of the overall cosmological Dark 
Matter. 

4.2 M2 dependence 
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Figure 4: (a) AiV 12 , as given by eq.(jl]) and (b) the relative size of radiative corrections 
to the relic density of Xi: as defined in eq.((7]), against M 2 , for our benchmark scenario 
© with Mi = 400 GeV. 

In Fig. Hk, b we plot the dependence of the absolute correction to the wino coupling 
of Xi an d the relative correction to the Xi renc density, respectively, on the wino mass 
parameter M 2 for bino mass parameter Mi = 400 GeV. Note that N i2 < here, i.e. the 
corrections increase the wino coupling |iVi 2 |. The bulk of these corrections is RG-like, i.e. 
can be attributed to the running of the SU (2) component of the xf f "g au g e " coupling. 
At first glance it might be surprising that the size of these 0(g 2 ) corrections is somewhat 
smaller than the O(gy) correction shown in Fig. |3k, even though g 2 ~ 3.3^. The 

§The higgsino components of these gaugino-like states, and the gaugino components of the higgsino- 
like states, still differ between these two tree-level calculations, but these have little effect on the relic 
density in this region of parameter space. Presumably these differences become very small if one compares 
one-loop calculations, since the two "tree-level" neutralino mass matrices only differ by 0(a) terms. 
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reason for this is that in the on-shell scheme the wino coupling receives a sizable negative 
"constant" correction (which does not grow logarithmically with the sfermion masses): 
the rather small value of sin 2 9w — 0.222 leads to a tree-level £77(2) gauge coupling 
g — e/ sin#iy that is somewhat larger than the coupling used in standard micrOMEGAs, 
where sin 2 6 W = 0.231. The total wino coupling g(iV 12 + AiV 12 ) is independent of the 
choice of renormalization scheme. 

Results for AUn and AVu, which describe the wino couplings of the lighter chargino, 
are very similar in magnitude to AiVi 2 ; these chargino couplings contribute to certain 
co-annihilation reactions. 

At M 2 ~ 390 GeV there is a discontinuity in AiVi 2 . This is because the dominantly 
wino-like neutralino no longer remains the lightest one, the bino-like one becomes the 
lightest state; we had found an analogous discontinuity in Fig. [3J As M 2 is increased 
beyond this value, | A7V 12 1 quickly decreases, simply because \N 12 \ decreases as Xi becomes 
increasingly bino-like. 

In most of the region with wino-like Xi that satisfies the LEP constraint M 2 > 110 
GeV [46] the corrections to the Xi re hc density shown in Fig. |4b are quite small and 
positive. Here both Xi annihilation into W + W~ pairs, and various Xixt co-annihilation 
reactions are important. As in Fig. [3b for Mi > 400 GeV, co-annihilation into leptonic 
final states is strongly suppressed by destructive interference between W ± and slepton 
exchange diagrams. The increase of the magnitude x? an d Xi couplings to £77(2) dou- 
blet (s)leptons further strengthens this cancellation; since this final state is in any case 
subdominant, the corrections increase the predicted relic density only slightly. 

Within the region where Xi is wino-like, the corrections to the relic density become 
significant only in the vicinity of the s— channel A, H poles at M 2 ~ 250 GeV. The cou- 
plings of these Higgs bosons are essentially proportional to iV 12 here, N\\ being very small. 
Since our corrections increase |iV 12 |, as noted above, they increase the Higgs exchange 
contributions to the annihilation cross section, and thereby reduce the relic density. Re- 
call also that our tree-level £77(2) gauge coupling is about 1.5% larger than that used 
by micrOMEGAs; this also contributes to the corrections as defined in eq.([7]). Note that 
there is also an H ± exchange pole in XiXi co-annihilation, which also contributes in this 
range of M 2 since the masses of A, H and H ± are within a few hundred MeV of each 
other. Since, as noted above, the wino coupling of xf is enhanced similarly as that of 
Xi, our corrections also enhance these contributions to the co-annihilation cross section. 
However, in both cases the Higgs couplings to x currents are suppressed by the small 
higgsino components of Xi andY]* 1 ; hence Higgs exchange in the s— channel is important 
only for a narrow range of M 2 |j 

Finally, for M 2 > 400 GeV x,i becomes increasingly bino-like, as mentioned above. 
Here the corrections to the bino coupling, AV n , are of similar size as those shown in 
Fig. |3k for Mi < 400 GeV. In the transition region, co-annihilation with x 2 an d Xi 

^Since a wino-like isn't an input state in our on-shell renormalization of the x sector, its mass 
receives non-vanishing one-loop corrections in our calculation. However, they amount to less than 200 
MeV in magnitude, significantly smaller than the total decay width of the heavy Higgs bosons, which 
amounts to ~ 1 GeV in our benchmark scenario. These corrections to the x? mass are therefore essentially 
irrelevant near the heavy Higgs poles, but lead to sizable changes of the relic density near the h pole; 
however, the lower bound on the \i mass from LEP experiments |46j excludes the region of parameter 
space where a wino-like Xi nas mass near rrih/2. 
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is important. For yet larger values of M2, near the end of the range shown in Fig. HI 
co-annihilation with these wino-like states is suppressed by the large bino-wino mass 
splitting, and annihilation into hA, ZH and leptonic final states dominates. Since our 
corrections increase the bino coupling of Xi by almost 2%, they increase many of the 
relevant cross sections, thereby reducing the predicted relic density by several percent. 

In standard cosmology, the Xi re hc density reaches its desired value for M 2 ~ 425 
GeV, where x^l still has significant wino-componentd, and co-annihilation with x° an d 
xf is still significant. Our corrections reduce the relic density by about 2% at this point. 
For M2 < 400 GeV, Xi could only contribute less than 10% of the desired Dark Matter 
density. 



4.3 fi dependence 

Having discussed our corrections for bino- and wino-like Xv we now consider a region 
of parameter space where x? is higgsino-like, by fixing M\ = 400 GeV, M 2 = 600 GeV 
and varying fi, chosing fi > 0. Note that there are two higgsino-like neutralinos. One of 
these states shows significantly stronger higgsino-gaugino mixing, which simultaneously 
reduces its higgsino content and its mass. Since the other, heavier higgsino-like state also 
has (even) larger higgsino content, its mass is one of the input masses in our version of 
the on-shell renormalization of the x sector. This means that the Xi mass does receive 
radiative corrections in our scheme, which can amount to several GeV. At the same time, 
the gaugino couplings of Xi receive corrections whose relative size is similar to those 
shown in figs. [3^, and H^i for the bino- and wino-coupling, respectively; recall, however, 
that the bino- and wino-couplings of higgsino-like states also receive corrections that can 
be interpreted as corrections to the off-diagonal entries of the x mass matrices, rather 
than corrections that come from the running of the "gauge" couplings of the x states. 
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Figure 5: Relative size of the corrections to (a) the mass and gaugino couplings of x? ; 
and (b) to the predicted Xi relic density against fi, for our benchmark point (|HJ) with 
Ah = 400 GeV and M 2 = 600 GeV. 



"Of approximately equal size, but opposite sign, of a photino state. 
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In Fig. [5^ we plot the relative size of the corrections to the mass and to the gaugino 
entries Nu and N12. For \i < 300 GeV the absolute correction to the Xx mass is approxi- 
mately 1.3 GeV, nearly independent of /j,, leading to a falling relative size of the correction 
with increasing fi. For even larger values of [i the bino-component of \\ begins to become 
sizable; this leads to a rapid decrease of the correction to the mass, which turns negative 
at fi ~ 370 GeV. Note, however, that the correction remains small in magnitude even at 
/1 = Mi = 400 GeV, where other versions of the on-shell scheme are badly behaved [31] . 
For yet larger values of fi, Xi is bino-like; its mass is then an input, and does not receive 
any corrections. 

Much of the correction to the relic density shown in Fig. [5b can be explained through 
this correction to the x? mass. Note that higgsino-like neutralinos strongly (co-) annihilate 
into final states containing two gauge bosons^] co-annihilation into // final states via 
s— channel exchange of a gauge boson is also important, while t— and u— channel sfermion 
exchange diagrams are suppressed by large sfermion masses (in case of squarks) and/or 
small Yukawa couplings of the corresponding fermions (in case of first and second gener- 
ation sfermions). As a result, for most of parameter space none of these cross sections is 
affected significantly by our corrections to Xi couplings^ 

The total correction to the relic density can nevertheless be sizable. This is most 
obvious in the vicinity of the heavy Higgs poles, where the positive shift of the Xi mass 
is important, since it slightly exceeds the width of the heavy Higgs bosons. It therefore 
significantly increases the cross section below the resonance, where the correction reduces 
the mass gap to the resonance, but reduces it above the resonance, where the corrections 
move Xi further away from the resonance. The positive correction to the Xi mass also 
increases the annihilation cross section for rather light Xi> close to the W + W~ and 
ZZ thresholds where the corresponding cross sections are still significantly phase space 
suppressed. On the other hand, away from poles and thresholds increasing the Xi mass 
slightly decreases the annihilation cross section, which basically scales like l/ m |o m this 
region. 

In principle one can chose a renormalization scheme where the LSP mass is used as 
input in the on-shell renormalization of the x sector. However, such a scheme would 
lead to quite large corrections in some regions of parameter space [34]. Moreover, co- 
annihilation reactions play a significant role in determining the relic density of a higgsino- 
like LSP. If the LSP mass is kept fixed, there will be significant loop corrections to the 
masses of the other higgsino-like states, m^o and m~±J3 as a result, there would be large 
corrections to co-annihilation cross sections proceeding through exchange of a heavy 
(neutral or charged) Higgs boson. The total size of the correction to the relic density 

*In principle W + W~ and ZZ final states can be produced via the s— channel exchange of the CP-even 
Higgs bosons, whose couplings to two neutralinos or two charginos do receive some corrections. However, 
these contributions are suppressed by the small size of the gaugino components of the (co-) annihilating 
states; moreover, since m\ M§, the couplings of the heavy CP-even Higgs boson H to two gauge 
bosons is very small. In practice our corrections to Xi couplings therefore do not significantly modify 
the (co-) annihilation cross sections into two gauge bosons. 

^The corrections to the "top Yukawa" coupling of Xi> which is oc A./V14, does play a role in the tiny 
kink visible at the ti threshold, /j, ~ 180 GeV. 

•fA renormalization scheme were the masses of all three higgsino-like states were used as inputs would 
be very badly behaved indeed, since then no input mass would have significant dependence on Mi or 
M2, leading to very large counterterms for these parameters [34] . 
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near the heavy Higgs poles would then still be of similar size as in our scheme. 

In the vicinity of the Higgs poles, our corrections to the wino- and bino-couplings 
of Xi (and, in co-annihilation reactions, of as we U as the wino couplings of xf) do 
play a role. Note that the wino- and bino-components of Xi have opposite sign in this 
region of parameter space, which means that they add constructively to the XiXi{H, A) 
couplings. We see in Fig. that in the region of the heavy Higgs poles, the relative 
corrections to Nu and Ni 2 are both positive, leading to an increase of the magnitudes of 
the XiXi{H, A) couplings. For larger \i the correction to iVi2 turns negative, helping to 
reduce the overall size of the corrections to the relic density in that region of parameter 
space. 

The corrections to (the gaugino part of) Higgs couplings to x currents also play a role 
in the (co-) annihilation cross sections into one massive gauge boson and one heavy Higgs 
boson, which open around \x = 310 GeV, as well into the light Higgs boson h and one of 
the heavy Higgs bosons A, H, H^, which become accessible for /i > 322 GeV; the former 
mostly proceed via the exchange of a heavy Higgs boson in the s— channel, while the latter 
also receive significant contributions from the exchange a gaugino-like x states in the t— 
or u— channel. With the exception of the region with strongest bino-higgsino mixing, 
\i ~ 400 GeV, where the corrections to iV 12 are relatively large and negative, all relevant 
Higgs couplings are enhanced by our correction, thereby increasing the corresponding 
cross sections, whereas the increase of the mass still tends to reduce the cross section due 
to the proximity of the H, A poles. These two effects tend to cancel, reducing the overall 
size of the correction to the relic density to < 1.5% for 300 GeV < fi < 400 GeV. 

For yet larger /i, Xx becomes bino-like, as noted above. The absolute size of the 
correction to Q^o then increases with increasing /i, as co-annihilation reactions, which are 
little affected by our corrections, become less important. Since the corrections increase the 
bino coupling of Xi, they reduce the relic density. The correction is somewhat smaller than 
in most of Fig. [3U since even at fi = 600 GeV final states containing one massive gauge 
boson and one Higgs boson, as well as diagrams with Higgs exchange in the s— channel, 
play significant roles; in both cases only one vertex gets corrected. Note finally that the 
correction is perfectly well-behaved at fi = M 2 = 600 GeV, as advertised (and in contrast 
to the scheme used in [52]). 

Within standard cosmology, the x? re hc density reaches its desired value at \i ~ 460 
GeV, in the region of strong higgsino-bino mixing, where corrections to Q^o are negative. 
For smaller values of n, x\ would only contribute a fraction of the total Dark Matter 
density. 

4.4 Co- annihilation with f\ 

If the mass difference between the LSP and the next-lightest superparticle (NLSP) is 
sufficiently small, co-annihilation processes play a role in the determination of the LSP 
relic density [53]. This is generic for wino- as well as higgsino-like LSP, where Xi 1S 
automatically close in mass to Xx'i in case of higgsino-like LSP, x 2 * s a ^ so nearby. In this 
subsection we study co-annihilation of a bino-like LSP with the lightest stau state f\. 
This is well motivated even within the framework of the minimal supergravity (mSUGRA) 
model [51]. The effect of co-annihilation can be described quite well [53] by introducing 
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an effective annihilation cross section, which is essentially the sum of the usual LSP 
(self-) annihilation cross section and the Xi^i and f\f\ annihilation cross sections, where 
the latter are suppressed by one and two powers, respectively, of the Boltzmann factor 
exp[— Tp/(m fl — m%o)]. Since the freeze-out temperature T F is basically proportional to 
m^o, the Boltzmann suppression factor depends on the relative mass splitting between 
the LSP and the lighter stau. 

We have assumed \x = 600 GeV and M 2 = 400 GeV. All other input parameters 
are as described in ([8]), except that we have used rriA = 1.0 TeVE, rn? R = Mi + 5 
GeV and tan (3 = 5. Note that the physical f\ mass is increased by the "D— term" 
contribution — sin 2 6 W cos(2/3)M§ > 0, and reduced by f L — f R mixing. For the given 
choice of parameters, we find a fi — Xi mass splitting of about 14 GeV, almost independent 
of Mi in the range shown in Fig. |HJ 




Figure 6: (a) Correction to the Xi annihilation cross section and (b) radiative correction 
to the relic density of Xi (i n %) against M\. Results are for M 2 = 0.4 TeV, \i = 0.6 TeV, 
rriA = 1.0 TeV, tan/3 = 5 and m fR = Mi + 5 GeV; the remaining parameters are as in 
our benchmark scenario (151) . 

In Fig. O (a), the correction to the Xi annihilation cross section is shown. Since fi is 
much lighter than all other sfermions in the scenario considered here, Xi pairs dominantly 
annihilate into t + t~ pairs via fi exchange in the t— or u— channel. Both vertices in these 
diagrams receive significant positive corrections, leading to an enhanced Xi annihilation 
cross section. 

*We have increased tha in order to avoid the heavy Higgs poles, where micrOMEGAs (version 2.2.CPC.i) 
produces a spurious divergence in the co-annihilation cross section. This occurs because f\ — > r + xi 
can occur on shell. As a result, fi + \i —> T + {A,H) can proceed through "exchange" of an on-shell 
Xi, giving an infinity in the (differential) cross section; note that Xi is stable, and hence has vanishing 
width. This divergence is not physical. In the present context, the very fast reactions fi O Xi + T help 
to maintain the relative equilibrium between the Xi an< l densities at temperatures well below Tp. 
In the formalism of ref. [53] these fast reactions are integrated out to arrive at the effective Boltzmann 
equation for the Xi density. The removal of this unphysical divergence is a separate issue not related 
to the radiative corrections we are interested in, so we simply avoid it by increasing the masses of the 
heavy Higgs bosons. 



24 



At M\ ~ 175 GeV the relative size of the correction begins to diminish somewhat. 
At this point annihilation to a tt pair opens. Increasing M\ along with rrif R reduces 
the annihilation cross section into t + t - pairs roughly like the inverse square of the Xi 
mass, but increases several other cross sections. In particular, since fi is kept fixed, 
increasing Mi increases higgsino-gaugino mixing, which leads to enhanced cross sections 
for W + W~ , ZZ as well as Zh final states. In some of these channels, as well as for the tt 
channel, there is destructive interference between s— channel Higgs exchange and t— and 
u — channel diagrams. The corrections to the annihilation cross sections into W + W~ and 
tt final states are therefore negative (about —0.5% and —3%, respectively). This leads to 
a reduction of the relative size of the correction to the total Xi annihilation cross section. 

The right frame of Fig. [HI shows that the correction to the x? ren c density is signifi- 
cantly smaller in magnitude than that to the Xi annihilation cross section. This is mostly 
due to co-annihilation contributions. Recall that the importance of these contributions 
depends on the relative f\ — x\ mass splitting; since the absolute mass splitting remains 
approximately constant over the entire range of M\ shown, the relative mass splitting 
decreases approximately oc 1/M\. As a result, co-annihilation becomes more important 
at larger M\. 

The most important co-annihilation reactions are X^i^t ~~ * T± + (7? where only 
one vertex is corrected, as well as f^f-f —> 77, jZ, which proceed solely through (true) 
gauge interactions, which we do not correct. However, f^f^ — > r ± r ± also contributes 
significantly, where again both vertices of the dominant x° exchange diagrams get cor- 
rected. Nevertheless, our corrections affect the co-annihilation cross sections much less 
than the x? annihilation cross section. Together with the increased importance of these 
co-annihilation contributions this explains the decline of the relative size of the correc- 
tions with increasing M\. In the case at hand, co-annihilation is dominant for M\ > 150 
GeV. For much larger values of M\ , the Boltzmann suppression factors are close to unity, 
and the relative size of the corrections stays approximately constant until one approaches 
the heavy Higgs poles (not shown). 

Within standard cosmology, in the parameter space shown in Fig. El the x\ relic 
density is three to four times higher than the desired value. This could be corrected by 
also reducing the masses of some other sleptons, and/or further reducing m ?1 . In both 
cases the importance of co-annihilation contributions would increase, so that the relative 
corrections to the relic density would be similar to that in the high— Mi region of Fig. [61 

5 Conclusions 

In this paper we have studied the impact of the effective couplings, as introduced in [32] . 
and of corrections to the neutralino and chargino masses, in the scheme introduced in 
ref . [31] , on the calculation of the cosmological relic density of the lightest neutralino in 
the MSSM. 

In Sec. 2 we analyzed the nature of the corrections captured by these effective cou- 
plings. We pointed out that they not only describe the difference in the running of true 
gauge couplings and gaugino couplings, but also corrections to the entries of the chargino 
and neutralino mass matrices that mix higgsino and gaugino states. Technically, these ef- 
fective couplings include x two-point function corrections plus associated counterterms. 
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These corrections are thus enhanced by multiplicity factors, since all (s)fermions con- 
tribute similarly to the gaugino two-point functions. It has also been known for some 
time that the magnitude of these corrections increases logarithmically with the mass of 
heavy sfermions [25] • There is thus some reason to believe that they capture a large 
fraction of all electroweak radiative corrections; in case of x° ~~ ffXi decays this has 
been shown to be the case [33] . 

In Sec. 3 we briefly described how we incorporated these corrections into the micrOMEGAs 
program package [36] . The corrections apply to both the "gauge" and the "Yukawa" com- 
ponents of xff couplings, as well as to the "gaugino" components of Higgs couplings to 
X currents. It should be noted that we use on-shell renormalization in the electroweak 
sector, whereas micrOMEGAs uses MS inputs for the QED coupling and the weak mixing 
angle. We nevertheless compare with micrOMEGAs tree-level predictions when illustrating 
the effect of the corrections, see eq. 

Numerical results are shown in Sec. 4. We found that the corrections to o~(x1x1 ~> 
e + e~) can reach +15% for wino-like xi- However, in this case Xi dominantly annihilates 
into W + W~ final states, which are relatively little affected by our corrections (only via 
s— channel Higgs exchange diagrams). The corrections to the total xi annihilation cross 
section due to the effective couplings therefore never exceed 4% for our choice of squark 
masses near 1.5 TeV, i.e. not far beyond the current LHC lower bound. The correction 
to the x^i mass, which never exceed a couple of percent in our scheme, can lead to 
significantly larger corrections to the total annihilation cross section in the vicinity of 
Higgs poles. These mass corrections were shown to be especially important for higgsino- 
like LSP. In all cases the relative corrections to the relic density are slightly smaller than 
those to the (effective) annihilation cross section, since an increase of the latter also 
reduces the temperature at which xi decouples from the thermal plasma. 

It should be noted that the corrections to the relic density are frequently comparable 
to, or larger than, the present observational uncertainty on the total Dark Matter density 
(in the framework of standard cosmology); this uncertainty is expected to be reduced 
significantly as soon as the PLANCK collaboration completes the analysis of their data 
on the cosmic microwave anisotropies. 

Note also that the effective couplings as defined in [32] include certain diagrams con- 
tributing to a complete one-loop calculation of o-(x1x1 ~~ >* f f) exactly, namely the two- 
point function corrections on the external xi lines with fermion-sfermion loops. Their 
inclusion into an "effective tree-level" calculation, as done here, can therefore be under- 
stood as a first step towards a full electroweak one-loop calculation [16] , which is however 
computationally far more complex. We use the same effective couplings also when a x is 
exchanged in the t— or u— channel, which contributes to sfermion co-annihilation; since 
here x is off-shell, the effective couplings include two-point function corrections to its 
propagator only approximately, but the relation with the exact calculation is still fairly 
straightforward. 

The situation is a bit more complicated for the Higgs couplings to x currents. Here 

*The corrected couplings g(N + AN), g{U + SU), g(V + AV) are scheme-independent up to two- 
loop terms; however, the tree-level couplings gN etc. are scheme-dependent already at one-loop order. 
We use micrOMEGAs choices here since this is currently the industry standard for LSP relic density 
calculations. 
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we correct only the gaugino part of the coupling. In the limit of vanishing gaugino- 
higgsino mixing, this corresponds to including only the gaugino to gaugino and higgsino to 
gaugino two-point functions, but no higgsino to higgsino or gaugino to higgsino diagrams. 
However, in realistic situations with non-vanishing \ mixing this no longer corresponds 
to a complete subclass of diagrams. Moreover, it has been shown some time ago [S5] 
that at least some corrections to ZxiXi couplings can be described in terms of neutralino 
two-point function corrections. There might thus be further scope for extension of the 
effective coupling scheme to Higgs couplings, and/or to couplings to W and Z bosons to 
X states. We plan to investigate these issues in future. 
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A Weak Sector 

We apply on-shell renormalization of the electroweak sector, where Mw and Mz are 
physical (pole) masses, and cos 9w = Mw/Mz- This gives [56] : 

5M W = mtE ww (M^); 

5M 2 Z = mtZzziMz); 

. a M w fSM w 5M Z \ fn . 

dcos9 w = — — — — — . (9) 

M z \M W M z J K J 

Here T>ww an d £zz are the transverse components of the diagonal W and Z two-point 
functions in momentum space, respectively; only the real parts of the loop functions 
should be included as indicated by the £Hc symbols. Note that, while the definitions of 
these three counterterms are formally as in the SM, in the MSSM there are many new 
contributions to Y* ww and Y, Z z involving loops of superparticles and additional Higgs 
bosons. 

In our treatment we have to calculate the counterterms for Mw, Mz and cos6w twice: 
for the purpose of computing effective couplings, we only include loops of matter fermions 
and sfermions, but for the correction to neutralino and chargino masses [34] we include 
the full MSSM contributions. 

The SU(2) gauge coupling is given by g = e/ sin 9 W , where the e is the QED couplings 
constants, the renormalization of which is discussed in the following Appendix. The 
counterterm 5g is then given by 

. 5e SsinOw Se cos 9w . . 

°9 = — a 9^~ a — = ■ a +9 . 2 n ocos9 w , 10 

sm t>w sin &w sin t) w sin 9 W 
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where 5cos6\y has been given in eq.([9]), and 5e = e(Mz)5Z^ Mz ^ with bZ e } Mzl given in 
eq. fl2"U|) below. The counterterm 5g is needed only in the calculation of effective couplings, 
so only matter (s)fermion loops are included. 

Similarly, the counterterm 5 tan 6w needed in eqs.(j4]) is to one-loop order given by 

x + a 5cos6 w 

dt&n9 w = — — — — . (11) 

sm V w cos Uw 

B Charge Renormalization 

The on-shell renormalization condition for the electric charge requires the electric charge 
to be equal to the full eey coupling for on-shell external particles in the Thomson limit. 
This leads to 

e barc ^ e ( Q ) p + 6z (0)j _ ( 12 ) 

The renormalization constant is adjusted to cancel the loop corrections to the eej vertex 
in this limit. This gives [57] : 

5Z ^ = ^ E -r7(? 2 )l« 2 =o + tan^^P, (13) 

where the £ again refer to the transverse part of the corresponding two-point func- 
tions. We can thus identify the renormalized charge e(0) = a/47to(0) where a(0) = 
1/137.036 ... is the fine structure constant, as defined in the Thomson limit. 

While pure on-shell renormalization of e is perfectly fine in a theory of photons and 
electrons, a problem arises in the presence of light quarks. The counterterm 5Zf^ receives 
contributions from loops involving a light fermion / that scale oc log(m^), leading to terms 
proportional to a log (m^/s) in corrections to physical processes at some energy scale s; 
these can be understood as describing the running of the QED coupling. The problem 
is that the masses of the light u, d, s (and, to some extent, c) quarks aren't well known, 
and in fact are difficult to define in the Thomson limit. 

Fortunately this dependence on the light quark masses can be traded for a dependence 
on experimentally measured cross sections for e + e~ annihilation into hadronic final states. 
To that end, one writes 

flight / in loops (g2) | g2=Q = Aa+ -L^eEijf / - 1°°PS(M|) , (14) 

where Aa is a finite quantity; moreover, the second term in eq.( IT4|) is finite for rrif — > 0, 
so that the u, d, s masses can be neglected when evaluating this term. Aa can be split 
into two parts: the contribution from the e,/x, r leptons and the contribution from the 
light quarks (i.e. all except t), Aa = Aa\ cpt + Aa had . Aa lcp has been calculated up to 
3-loop order as [58] 

Acti cp = 0.031497687, (15) 
while Acthad can be extracted from experiment via a dispersion relation [59] : 

Aa had = 0.027626 ± 0.00138. (16) 
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This leads us to 

SZ^ = lA^allloops^l^^^ght/inloops^l^J 

+ l Aa + f m l0 ° PS(M|) + tan % ^MT" (18) 

One can achieve even higher accuracy by re-summing the light fermion correction. 
This corresponds to using the "running on-shell" coupling a(M§) = ct(0)/(l — Aa) = 
1/128.93 as input ("tree-level") coupling[j] The corresponding renormalization constant 
becomes: 

e bare ^e(0)(l + ^°)) = e(0) (l + ^ Aa - ~Aa) (l + *Z?>) 

= e(M|) (l + 6Ze (M2z) + higher order) . (19) 

Thus, 

2 

+ 2^|^ ht /inlOOPS ( M i)+ tan ^^T- ( 2 °) 

We need 5Zt^ Mz ^ only in the calculation of effective couplings, hence only the matter 
(s)fermion loop contributions are included here. 



C Higgs Sector 

The only relevant counterterm determined from from this sector is 5tan/3. It is fixed by 
forbidding transitions between the CP-odd Higgs boson A and the Z boson on the A 
mass shell [601 ED E2]. This gives: 

* tan0 =2MOT 3m( ^ (mi,) - <21) 

Note that the couplings of A contain an extra factor of imaginary % relative to the gauge 
couplings. Therefore the imaginary part of the two-point function appears in eq. (l2~TT) : 
this contains the real (dispersive), infinite, part of the loop function. 

The counterterms for cos {3 and sin f3 are thus: 

<5cos/3 = -cos 3 /3tan/35tan/3, (22) 

^As mentioned in Sec. 3, this differs slightly from the running MS coupling a(Mz) = 1/127.944 often 
used in tree-level calculations of electroweak processes; see e.g. the review by Erler and Langacker in 

US- 
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8 sin (3 = cos 3 (38 tan (3. 



(23) 



The fermion loop contribution to E^z is proportional to the mass of the fermion in the 
loop. In case of SM quarks and leptons, this gives an additional factor of the corresponding 
Yukawa coupling times M w /g. These contributions to 8 tan (3 should therefore be counted 
as 0(X 2 ), not as O(Xg), the second factor of A stemming from the coupling of A to the 
SM fermion in the loop. 

These counterterms are needed both for the calculation of one-loop corrected neu- 
tralino and chargino masses and for the evaluation of the effective couplings. As in case of 
the counterterms to the W and Z masses, in the former case all MSSM loop contributions 
are included, while in the latter case only matter (s)fermion loops are included. 



D Chargino and Neutralino Sector 

In this appendix, we describe the determination of the relevant counterterms in the 
chargino-neutralino sector of the MSSM, closely following the discussion in [32 | |35 | |63 | 134"] . 

The one-loop mass eigenstates can be related to the tree level mass eigenstates through 
wave function renormalization: 

X* bare = (Sij + \8Z^P L + \8ZfjPn) x y™ alized . (24) 

where Pl and Pr are the chiral projectors. This relation holds for both the chargino 
sector, with Xi = xti * ^ {1)2}, and the neutralino sector, with Xi = Xii i £ {1)2,3,4}. 
The entries of the mixing matrices N, U and V have been chosen to be identical to their 
tree level values [53]. 

The corresponding masses receive explicit corrections from one-loop diagrams, which 
can be written as 
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8m f = ~m f ^KeEj/(mJ) + <HeE}/(m})j + - ^KeEj, (mj) + 9teE^(m^)J , (25) 

where / denotes a fermion species with mass m/. 91c again denotes the real parts of the 
loop integrals involved and the E refer to various terms in the general two-point function 
connecting fermions / and /' in momentum space: 

Z ff ,(p) = f [P L ^{p) + P R ^ V fp{p)] + P L Z S f f,(p) + P R X s f $(p) ■ (26) 

Only the case / = /' is needed in eq. ([2"5"j) . but we will need off-diagonal two-point 
functions for the determination of the wave function counterterms. 

The corrections of eq. (]25|) are in general divergent. These divergencies are absorbed 
into counterterms to the mass matrices M c and M n of the charginos and the neutralinos 
respectively, i.e. 

M barc = M rcnormalizcd + 8M. (27) 

As usual, we define the physical (on-shell) masses as poles of the real parts of the (one- 
loop corrected) propagators. The physical chargino masses are then given by, 

m « + = m + + (U*8M c V- l ) u - 8m~ + ; (28) 
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the corresponding expression for the neutralinos is 

m ° s = m -o + (N*5M Q N- 1 ) ii - . (29) 

The masses m~+,o appearing on the right-hand sides of eqs. f[28]) and fl29|) are the (finite) 
tree-level masses. 5m~+,o are the explicit loop corrections of eq. ([25]) as applied to the 
charginos and neutralinos. Finally, 5M c,n are the counterterm matrices of eq. (I27p for the 
chargino and neutralino sector, which we yet have to determine. 

To one-loop order counterterms to products like My/ sin /3 can be written as (SMw) sin (3+ 
M\yS sin j3, and so on. 

Altogether, the tree-level chargino and neutralino mass matrices depend on seven 
parameters, hence there are seven different counterterms: 5M\y, 5Mz, 5@w, $ tan /3, 5Mi, 
5M 2 and 5fi. The first three of these already appear in the SM, see Appendix A. <5tan/3 
has been fixed in the Higgs sector, see Appendix O 

Hence only the counterterms to M±, M 2 and \i remain to be fixed in the chargino 
and neutralino sector. This means that we can only chose three of the six masses in 
this sector to be input masses which are not changed by loop corrections, i.e. only three 
of the six "tree-level" masses are physical (all-order) masses. We choose the wino-like 
chargino mass, the bino-like and higgsino-like neutralino masses as inputs (i.e. physical 
masses). Ref.JM] explains in detail why this particular choice of input states leads to the 
most stable perturbative expansion, by avoiding spurious large corrections. The other 
three masses will receive finite, but non-zero corrections. The remaining counterterms 
can then be computed from eqs. (l2"gj) and (T2"9"|) : 



SMi 
5M 2 
5fi 



D 

Nm 2 
D ' 

D ' 



(30) 
(31) 
(32) 



where 

D = 2N; 3 N* 4 Nf 1 U* kl V: i 



2iv* N* 3 N;,u* kl v k \ + n; 2 n*, u* k2 v k * 2 - iv* n; 2 u* k2 v k * 2 , (33) 



N„. 



N*, ( Smg - 25M^ 3 N* l N* 3 - 26M^ 3 N* 2 N* 3 - 28M^N* l N* i - 26M^N* 2 N* 



-N*l (-Sm^o + 25Mf 3 iv; 1 iv; 3 + 25M^N* 2 N; 3 



25M^N* 1 N* i + 25M^N; 2 N; / , 



+ (-NgN* + NgNg) (-5m~ t + 5M^U* k2 V k \ + 5M^ 2 U* kl V k * 2 ) , 



u, 



kl 



N. 



M 2 



+ 

+ 
+ 



2^+^^^ - 25m- t N* 1 N* 3 N* 4 - 25 M^N^N*, U* k2 V k \ 
25M° 1 N£N* 3 N* 4 U* k2 V k \ - 26M c l2 N* 3 N* 4 Nf 1 U* kl V k * 2 - SmqNgUtffe 
SM c 12 N*;N; 3 N; A U* k iV k * 2 + Sm^N;tU* k2 V: 2 - 26M? 3 N* l N* 3 Nf 1 U* k2 V k * 2 
26M2 3 N* 2 N* 3 Nf l U* k2 V k * 2 - 25M? i N? 1 N* i Nf 1 U* k2 V k * 2 - 25M^N; 2 N;,Nf 1 U* k2 V k * 2 
26M? 3 N*;N* l N; 3 U* k2 V k * 2 + 25M^ N* 2 N* 3 U* k2 V k * 2 + 28M? A Nf l N* l N* A U* k2 V: 2 
26M^N* 2 N* 4 U* k2 V k * 2 , (35) 
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N 



- 25m- t N: 3 N* 4 Nf 2 - 25m^+ N* 2 N* 3 N* A - 25m^N* 3 N*, Utffo - 8m^N* 2 U* k2 V k \ 

- 25M? A N* 1 N* A N* 2 U* k2 V k * 2 - 28M^N* 2 N* A Nf 2 U* k2 V k % + 25M^N* 2 N* l N*JJi 2 V k \ 
+ 25M^N* 2 N; 2 N* 3 U* k2 V k * 2 + 25M? A N* 2 N* 1 N* A U* k2 V k \ + 25M% A N* 2 N* 2 N* A U* k2 V k * 2 

- 25M c 12 N* 3 N* A Nf 2 U* kl V k * 2 + 25M c l2 N* 2 N* 3 N* A U* kl V k * 2 + 45 M? 3 N* 3 N* A NJ.NJ^V^ 

+ 45MZ 3 N* 3 N*,N* 2 N* 3 U* kl V k \ - 25M? 3 N* 1 N* 3 Nf 2 U* k2 V k * 2 - 25M^ 3 N* 2 N* 3 Nf 2 U* k2 V k * 2 
+ 4 5 M™ 4 N* 3 N* A N* x N* A U kl V k \ + 4 6 M 2A N* 3 N* A N* 2 N* A U kl V kl + 25m^N* 3 N*M l V; i 

46M2 A N* 2 N* A N* 3 N* A U* kl V k \ - 25M c 21 N* 3 N* A Nf 2 Ut 2 V k \ + 25M c 2l N* 2 N* 3 N* A Ut 2 V k \ 
45M? 3 N* x N* 3 N* 3 N* A U* k iV k \ ~ ^M 23 N* 2 N* 3 N* 3 N* 4 U kl V kl - 4<5M 1 r 4 iV* 1 N* A N* 3 N* A U kl V kl 
+ 5m^Nf 2 U* k2 V k * 2 . (36) 

Here i and j are the indices of the two neutralino input states, and k is the index of 
the chargino input state; the counterterms are symmetric under % ■<->■ j, both D and 
the numerators being antisymmetric. Since we take input states with fixed physical 
characteristics (bino-, wino- or higgsino-like) , the numerical values of i, j, k will depend 
on (the ordering of) the parameters Mi, M2 and //; see ref.JM] for further details. 

These counterterms are needed both for the calculation of the one-loop corrected 
chargino and neutralino masses, where all MSSM loop contributions have been included 
in the evaluation of eq. (1251) . and for the calculation of effective couplings (see below), 
where only matter (s)fermion loop contributions have been included. 

The diagonal wave-function counterterms are determined by normalizing the residues 
of the propagators [351 E2] : 

6Zjf = -^^ VR {ml)-ml (^l L '(mi) +*e£™K)) -2m*KeI]f 

(37) 

The prime denotes derivative with respect to the external momentum squared. This ex- 
pression holds both for charginos [where these counterterms had an additional superscript 
+ in eqs.(j4])] and neutralinos. 

The off-diagonal wave-function renormalization counterterms are determined by de- 
manding that the corresponding off-diagonal two-point correlation functions vanish for 
on-shell external momenta [351 [32]. This givesQ 

L 2 (m^Ef(m|.) + m*.XfV|.) + m^m^E^(m|.) + m^Eg 1 ^)) 
5Z H = m 2 _ m 2 > ( 38 ) 

R _ 2 {m^tf{m\,) + m^tf(ml.) + m^m^.Eg L (m|.) + m^Eg^m?.)) 
5Z H = m 2 _ m 2 > ( 39 ) 

where, 

±l L > VR = &E™«™ f%> = *eSf - Srh tJ , Ef = £eEf - Srh*, (40) 



■^Note that our E,j(m|) corresponds to S Jl (m?) of FormCalc. 
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These expressions again hold for both charginos and neutralinos, where the latter were 
labeled with an additional superscript + in eqs.fll]) in the main text. Here Sfhij denotes 
(U*5M c V~ 1 )ij and (N*5M n N~ 1 )ij for charginos and neutralinos, respectively. The cal- 
culation of the counterterm matrices SM C and 5M n has been discussed in the first part 
of this Appendix. Hermiticity implies the following relations: 

sgV) = W) = s f V), W) = 4V) = 

Neutralinos are Majorana fermions. Consequently the following additional relations [31)1 
[63] hold: 

These relations can be used to further simplify eqs. (1401) as applied to the neutralinos 
[351 [63]. 

Note finally that neutralino couplings may be purely imaginary even if CP is con- 
served; this is true jbr states corresponding to negative eigenvalues of the neutralino 
mass matrix. Since fHc only acts on the loop functions, these imaginary parts have to be 
kept in eqs. fl371) - fl39l) : in fact, they are necessary to yield finite results in eqs. fl38l) . fl39|) 
for m^o = m^o, which is possible if M\M.2 < 0. In this case Xi an d can be com- 
bined into a Dirac neutralino state, and the corresponding "off-diagonal" wave function 
renormalization becomes part of the diagonal wave function renormalization of this Dirac 
neutralino. Degenerate charginos can only occur if tan/3 = 1, which is excluded by Higgs 
mass bounds, or in the presence of nontrivial CP-violating phases in the chargino mass 
matrix. 
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